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Abstract 

The tunneling of a giant spin at excited levels is studied theoretically in 
mesoscopic magnets with a magnetic field at an arbitrary angle in the easy 
plane. Different structures of the tunneling barriers can be generated by 
the magnetocrystalline anisotropy, the magnitude and the orientation of the 
field. By calculating the nonvacuum instanton solution explicitly, we obtain 
the tunnel splittings and the tunneling rates for different angle ranges of 
the external magnetic field {9h = it/2 and tt/2 < Oh < it). The temperature 
dependences of the decay rates are clearly shown for each case. It is found that 
the tunneling rate and the crossover temperature depend on the orientation 
of the external magnetic field. This feature can be tested with the use of 
existing experimental techniques. 
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I. INTRODUCTION 



Recently, nanospin systems have emerged as good candidates to display quantum phe- 
nomena at a mesoscopic or macroscopic scale.0 Theoretical investigations showed that quan- 
tum tunneling was possible in ferromagnetic (FM) nanoparticles containing as much as 
10 5 — 10 6 spins.B At extremely low temperature, the magnitude of the total magnetization 
M is frozen out and thereby its direction becomes the only dynamical variable. In the 
absence of an external magnetic field, the magnetocrystalline anisotropy can create ener- 
getically degenerate easy directions depending on the crystal symmetry. Tunneling between 
neighboring states removes the degeneracy of the original ground states and leads to a level 
splitting. This phenomenon is called macroscopic quantum coherence (MQC). However, 
MQC is hard to be observed in experiments without controlling the height and the width 
of the barrier. It has been believed that a magnetic field is a good external parameter to 
make the quantum tunneling observable. By applying a magnetic field in a proper direction, 
one of the two energetically equivalent orientations becomes metastable and the magne- 
tization vector can escape from the metastable state through the barrier to a stable one, 
which is called macroscopic quantum tunneling (MQT). A large number of experiments in- 
volving resonance measurements, magnetic relaxation, and hysteresis loop study, Mossbauer 
spectroscopy, and neutron scattering study for various systems showed either temperature- 
independent relaxation phenomena or a well-defined resonance depending exponentially on 
the number of total spins, which supported the idea of magnetic quantum tunneling.© 

To our knowledge, the tunneling of a single spin degree of freedom was first studied by 
Korenblit and Shender in 1978 JH More recently, the tunneling problem of the magnetization 
reversal was studied extensively for the single-domain FM nanoparticles in a magnetic field 
applied at an arbitrary angle. This problem was studied by Zaslavskii with the help of 
mapping the spin system onto a one- dimensional particle system.! For the same system, 
Miguel and Chudnovskyi calculated the tunneling rate by applying the imaginary-time path 
integral, and demonstrated that the angular and field dependences of the tunneling exponent 



1 



obtained by Zaslavskii's method and by the path-integral method coincide precisely. Kim 
and Hwang performed a calculation based on the instanton technique for FM particles 
with biaxial and tetragonal crystal symmetryjl and Kim extended the tunneling rate for 
biaxial crystal symmetry to a finite temperature.! The quantum-classical transition of the 
escape rate for uniaxial spin system in an arbitrarily directed field was investigated by 
Garanin, Hidalgo and Chudnovsky with the help of mapping onto a particle moving in a 
double-well potential.! The switching field measurement was carried out on single-domain 
FM nanoparticles of Barium ferrite (BaFeCoTiO) containing about 10 5 — 10 6 spins. i The 
measured angular dependance of the crossover temperature was found to be in excellent 
agreement with the theoretical prediction^ which strongly suggests the quantum tunneling 
of magnetization in the BaFeCoTiO nanoparticles. Lii et al. studied the quantum tunneling 
of the Neel vector in single-domain antiferromagnetic (AFM) nanoparticles with biaxial, 
tetragonal, and hexagonal crystal symmetry in an arbitrarily directed field.§ 

It is noted that the previous results of spin tunneling at excited levels in an arbitrarily 
directed field were obtained by numerically solving the equation of motion satisfied by the 
least trajectory! and the system considered in Ref. 6 had the simple biaxial crystal symme- 
try. The purpose of this paper is to present an analytical study of the quantum tunneling 
at excited levels in the FM particles with an arbitrarily directed field. Moreover, the system 
considered in this paper has a much more complex structure (i. e., the general structure in 
experiments), such as trigonal, tetragonal, and hexagonal crystal symmetry. By applying 
an arbitrarily directed magnetic field, the problem does not possess any symmetry and for 
that reason is more difficult mathematically. However, it is worth pursuing because of its 
significance for experiments and the easiest to implement in practice. Since the result of 
spin tunneling at excited levels for tetragonal symmetry is a generalization of that of tun- 
neling at ground-state levels studied by Kim and Hwang,! we can compare our results with 
theirs by taking the low-energy limit. We will show that MQC and MQT can be consec- 
utively observed by changing the direction of magnetic field, and discuss their dependence 
on the direction and the magnitude of field. The dependence of the crossover temperature 



T c and the magnetic viscosity (which is the inverse of WKB exponent at the quantum- 
tunneling-dominated regime T <C T c ) on the direction and the magnitude of the field, and 
the magnetic anisotropies is expected to be observed in future experiments on individual 
single-domain particles with different magnetocrystalline anisotropies. Both the nonvacuum 
(or thermal) instanton or bounce solution, the WKB exponents and the preexponential fac- 
tors are evaluated exactly for different angle ranges of the magnetic field (Oh = vt/2 and 
7r/2 < H < 7r). The low-energy limit of our results agrees well with that of ground-state 
spin tunneling. In order to compare theories with experiments, predictions of the crossover 
temperature corresponding to the transition from classical to quantum behavior and the 
temperature dependence of the decay rate are clearly shown in this paper. Both variables 
are expressed as a function of parameters which can be changed experimentally, such as 
the number of total spins, the effective anisotropy constants, the strength and orientation 
of applied magnetic field. Our results show that the distinct angular dependence, together 
with the dependence of the WKB tunneling rate on the strength of the external magnetic 
field, may provide an independent experimental test for the spin tunneling at excited levels 
in nanoscale magnets. When the effective magnetic anisotropy of the particle is known, our 
theoretical results give clear predictions with no fitting parameters. Therefore, quantum 
spin tunneling could be studied as a function of the effective magnetic anisotropy. Our 
results should be helpful for future experiments on spin tunneling in single-domain particles 
with different magnetocrystalline anisotropies. 

This paper is structured in the following way. In Sec. II, we review briefly some basic 
ideas of spin tunneling in FM particles. And we discuss the fundamentals concerning the 
computation of level splittings and tunneling rates of excited states in the double-well-like 
potential. In Sees. Ill, IV, and V, we study the spin tunneling at excited levels for FM 
particles with trigonal, tetragonal and hexagonal crystal symmetry in an external magnetic 
field applied in the ZX plane with a range of angles n/2 < H < 7r, respectively. The 
conclusions are presented in Sec. VI. 
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II. PHYSICAL MODEL OF SPIN TUNNELING IN FM PARTICLES 



For a spin tunneling problem, the tunnel splitting or the tunneling rate is determined by 
the imaginary-time transition amplitude from an initial state \i) to a final state |/) as 

U fi = (f\ e~ HT \i) = J VQexp (S E ) , (1) 

where S E is the Euclidean action and T>Vt is the measure of the path integral. In the 
spin-coherent-state representation, the Euclidean action is 



S E (6,((>) = j j dr 



\Mo(d<f>\ .Mq ( dy, , 
% — — — i — — I cosv + E 
7 \drj 7 \dr / 



(2) 



where V is the volume of the FM particle and 7 is the gyromagnetic ratio. Mo = |M| = 
hjS/V, where S is the total spin of FM particles. It is noted that the first two terms 
in Eq. (2) define the Berry phase or Wess-Zumino, Chern-Simons term which arises from 
the nonorthogonality of spin coherent states in the north-pole parametrization. The Wess- 
Zumino term has a simple topological interpretation. For a closed path, this term equals 
— iS times the area swept out on the unit sphere between the path and the north pole. The 
first term in Eq. (2) is a total imaginary-time derivative, which has no effect on the classical 
equations of motion, but it is crucial for the spin-parity effect s.0'0Hl3 However, for the closed 
instanton or bounce trajectory described in this paper (as shown in the following), this time 
derivative gives a zero contribution to the path integral, and therefore can be omitted. 

In the semiclassical limit, the dominant contribution to the transition amplitude comes 
from finite action solution (instanton or bounce) of the classical equation of motion. The 
instanton's contribution to the tunneling rate T or the tunnel splitting A is given by3 

T (or A)=Au; p (^-) 1/2 e- s «, (3) 



2ti t 

where uj p is the oscillation frequency in the well, S c i is the classical action, and the prefac- 
tor A originates from the quantum fluctuations about the classical path. It is noted that 
Eq. (3) is based on quantum tunneling at the level of ground state, and the temperature 
dependence of the tunneling rate (i.e., tunneling at excited levels) is not taken into account. 



However, the instanton technique is suitable only for the evaluation of the tunneling rate 
or the tunnel splitting at the vacuum level, since the usual (vacuum) instantons satisfy the 
vacuum boundary conditions. In this paper, we will calculate the nonvacuum instantons 
corresponding to quantum tunneling at excited levels. 

For a particle moving in a double- well-like potential U (x), the level splittings of degen- 
erate excited levels or the imaginary parts of the metastable levels at an energy E > are 
given by the following formula in the WKB approximation,^ 

AE (or Im E) = exp [-S {E)\ , (4) 

7T 

and the imaginary-time action is 



rx 2 (^) 

S (E) = 2v / 2m / dx^JU(x)-E, (5) 

Jx!(E) 

where X\^ (E) are the turning points for the particle oscillating inside the inverted potential 
—U (x). uj (E) = In jt (E) is the energy-dependent frequency, and t (E) is the period of the 
real-time oscillation in the potential well, 

t(E) = V2^r iE) / X (6) 

Jx 3 (E) y/E-U{x) 

where X3 4 (E) are the turning points for the particle oscillating inside the potential U (x). 



III. MQC AND MQT FOR TRIGONAL CRYSTAL SYMMETRY 

In this section, we study the quantum tunneling of the magnetization vector in single- 
domain FM nanoparticles with trigonal crystal symmetry. The external magnetic field is 
applied in the ZX plane, at an angle in the range of 7r/2 < 9h < 7r. Now the total energy 
E (6, 0) can be written as 

E {6, <f)) = K x sin 2 9 - K 2 sin 3 9 cos (30) - M H X sin 9 cos <p - M H Z cos 9 + E , (7) 

where K\ and K2 are the magnetic anisotropy constants satisfying K\ ^> K2 > 0, and Eq 
is a constant which makes E (9, 0) zero at the initial orientation. As the magnetic field is 



applied in the ZX plane, H x = Hs\nO H and H z = H cos Oh, where H is the magnitude of 
the field and Oh is the angle between the magnetic field and the z axis. 
By introducing the dimensionless parameters as 



K 2 = K 2 /2K 1 ,H X = H X /H ,H Z = H z /H , 



(8) 



Eq. (7) can be rewritten as 
1 



E 



sin — K 2 sin 9 cos (30) — H x sin 9 cos — H z cos 9 + E , 



(9) 



where E (9, 0) = 2K X E (9, 0), and H = 2K 1 /M . At finite magnetic field, the plane given 
by — is the easy plane, on which E (9, 0) reduces to 

1 



E(9,(f) = 0) = - sin 2 9 - K 2 sin 3 9 - H cos (9 - 9 H ) + E . 
2 



(10) 



We denote 9q to be the initial angle and 9 C the critical angle at which the energy barrier van- 
ishes when the external magnetic field is close to the critical value H c {6 H ) (to be calculated 
in the following). Then, the initial angle 9 satisfies [dE (0,4) = 0) /d0] e = 0, the critical 
angle C and the dimensionless critical field H c satisfy both [dE (0, = 0) /dOV 
and [d 2 E(0,4) = 0) /d0 2 ] ( 



~,H—H C 







,H—H C 



(sin^) 2/3 + |cos^| 2/3 
1 



0. After some algebra, H c (0 H ) and C are found to be 



1-3X2 



3/2 



(l + |cot^| 2 / 3 ) 1/2 



+QK 2 



(l + |cot^| 2/3 ) 



3/2 



sin 



1 + |cot^| 2/3 



2K 



I cot H \ 2/3 



|cot^-- |2/3 



7 H 



1/2 



4Kn 



I cot H \ 2/3 



|cot^- |2/3 



3/2 



(11a) 



• (Hb) 



Now we consider the limiting case that the external magnetic field is slightly lower than 
the critical field, i.e., e = 1 — H/H C <^ 1. At this practically interesting situation, the barrier 
height is low and the width is narrow, and therefore the tunneling rate in MQT or the tunnel 
splitting in MQC is large. Introducing r\ = C — (\r)\ <C 1 in the limit of e 1), expanding 
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[dE (9, = 0) /d9] e=6o = about 9 C , and using the relations [dE (9, = 0) /dO\ 0=e jj =1 j = 
and [d 2 E {9, = 0) /d9 2 ~\ Q _ Q jj =h - = 0, we obtain the approximation equation for r\ in the 
order of e 3 / 2 , 



-eH c sin (9 C -9 H )-rj 2 [- sin 29 c + 3K 2 cos 39 



+V 



eH c cos (9 C - 9 H ) + rj 2 [ - cos 29 c - 3K 2 sin 39, 



0. 



(12) 



Then E (9, 0) reduces to the following equation in the limit of small e, 



E (5, 0) = 2K 2 sin 2 (30/2) sin 3 (0 O + + sin (0„ + 5) (1 - cos 0) + £i (5) , (13) 
where 5 = 9 — 9 (\5\ <C 1 in the limit of e 1), and E 1 (5) is a function of only 5 given by 



Ei (5) = -- 
1 



H c sin (0 C — 9 H ) — K 2 ( cos 3 9 C — - sin 2 9 C cos 0< 



- [# c cos (9 C -9 H )- 3K 2 (sin 3 9 C - 4 sin0 c cos 2 C )] 
~K 2 (sin 3 c -4sin0 c cos 2 c ) 5 2 e. 



(35 2 V - 5 s ) 
3 



(14) 



In the following, we will investigate the tunneling behaviors of the magnetization vector 
at excited levels in FM particles with trigonal crystal symmetry at different angle ranges of 
the external magnetic field as 9 H = n/2 and n/2 < 9 H < ir, respectively. 

A. H = vr/ 2 



For 9 H = tt/2, we have 9 C = ir/2 from Eq. (lib) and rj = y^e (l + \K 2 ) from Eq. (12). 
Eqs (13) and (14) show that is very small for the full range of angles 7r/2 < 9n < it. 
Performing the Gaussian integration over 0, we can map the spin system onto a particle 
moving problem in the one-dimensional potential well. Now the imaginary-time transition 
amplitude Eqs. (1) and (2) becomes 



U 



fi 



dSexp (Se 



= Jd6e Xp \-Jdr \m{^) 



(15) 



with the effective mass 



hS 2 



m 



2VK 1 [1 - e + 9K 2 ] 



and the effective potential 



(16) 



The plot of the effective potential E\ (5) as a function of 5 (= 9 — 9 ) for 9 H = ir/2 is shown 
in Fig. 1, and HU (5) = 2K{VE\ (5). The problem is one of MQC, where the magnetization 
vector resonates coherently between the energetically degenerate easy directions at 8 — and 
5 = 2y / 2e (l + \K 2 ) separated by a classically impenetrable barrier at 5 = \[2e (l + \K 2 ). 

The nonvacuum (or thermal) instanton configuration 5 P which minimizes the Euclidean 
action in Eq. (16) satisfies the equation of motion 



(17) 



where E > is a constant of integration, which can be viewed as the classical energy of the 
pseudoparticle configuration. Then the kink-solution is 



$p = V + Vv 2 - asn (^i r > k) , 



(18) 



where a = 2^J j^y, and uj\ = ^^y/v 2 + a - sn(cjir, k) is the Jacobian elliptic sine function 



of modulus k = y^r^- The Euclidean action of the nonvacuum instanton configuration 
Eq. (18) over the domain (—(3,(3) is found to be 

2 







( dS p \ 


/ dr 













+ U{S P ) 



W + 2E(3, 



(19a) 



with 



W 



8 iWrnf 27 \ 1 
3V h V 2 7 ^1 - k' 2 /2 



E(k)- 



k 



12 



;K(k) 



(19b) 



2- k' 2 ' 

where k' 2 = 1 — k 2 , and (3 = 1/ksT with k% the Boltzmann constant. K (k) and E (k) are 



the complete elliptic integral of the first and second kind, respectively. The general formula 



Eq. (4) gives the tunnel splittings of excited levels as AE = ^p-exp(— W), where W is 
shown in Eq. (19b), and uj (E) = is the energy-dependent frequency. For this case, the 
period t (E) is found to be 

-5 2 



* (E) = [ 



d5 



2hm 



--K(k'), 



(20) 



l 5l y/E-U (5) "V^iV v /^T^ 

where 5\ = rj + \Jr\ 2 — a, and <5 2 = rj + ^/rf+~a. Now we discuss the low energy limit where 
E is much less than the barrier height. In this case, k' A = <C 1, so we can perform 



the expansions of K (k) and E (k) in Eq. (19b) to include terms like k' A and A;' 4 In (jj), 



E{k) 



1 

1 + 2 



lnl F 



ff<*)=ln(£)+i 



1 

2 



9 

k' 2 4- — 
k + 64 



k' 4 ■ ■ ■ . 



With the help of small oscillator approximation for energy near the bottom of the potential 



well, E n = (n + \) fix, Q ± = yJ±U" {5 = 0) = Eq. (19b) is expanded as 



1 - f - f if 2 / 1 
29/2^3/2 l n + 2 



(21a) 



where 



Wo = ^ (1 + 1 + 



(21b) 



Then the low-lying energy shift of n-th excited states for FM particles with trigonal crys- 
tal symmetry in the presence of an external magnetic field applied perpendicular to the 
anisotropy axis (0h = tt/2) is 

n+1/2 

exp (-Wb) . (22) 



hAE n = - _ 



For n = 0, the energy shift of the ground state is 

Ol3/4 / cr 7 \ 

£A£ = (i^) e^S^ 2 ( 1 - - + T ^ 2 J exp (-W ) . 
Then Eq. (22) can be written as 



(23) 



hAE n = % (hAE ) , (24a) 

where 

2 9/2 ^ 3/2 

* = 1 e 15^ • ( 24b ) 

Since we have obtained the tunnel splittings at excited levels, it is reasonable to study 
the temperature dependence of the tunneling rate. It is noted that Eqs. (24a) and (24b) 
are obtained under the condition that the levels in the two wells are degenerate. In more 
general cases, the transition amplitude between two levels separated by the barrier or the 
decay rate should be sensitive to this resonance condition for the two levels. If in the case of 
the potential with two degenerate levels only one of the levels is considered as a perturbative 
metastable state; however, a fictitious imaginary energy can be calculated by consideration 
of possible back and forth tunneling (i.e., by regarding the instanton-antiinstanton pair as a 
bounce-like configuration) in the barrier. Therefore there exists a relation between the level 
splitting and this imaginary part of metastable energy level, and has been referred to as the 
Bogomolny-Fateyev relation based on equilibrium thermodynamics^ 

lmE n = n (AE n ) 2 /Au (E n ) , (25) 

where u (E n ) is the frequency of oscillations at energy level E n . At finite temperature T the 
decay rate r = 2 Im E n can be easily found by averaging over the Boltzmann distribution 

r ( T ) = J" Im E « ex P (~ hE n(3) , (26) 

where Z = Xln ex P (~^nP) is the partition function with the harmonic oscillator approx- 
imated eigenvalues E n = (n + |) Oi. With the help of the Bogomolny-Fateyev relation Eq. 
(25), the final result of the tunneling rate at a finite temperature T is found to be 

r (T) = -L. (1 - e -"^) (AE ) 2 I (2 gi e-^) , (27) 
where AEq and q\ are shown in Eqs. (23) and (24b). I (x) = En=o ( x / 2 ) / ( 

nX) is the 

modified Bessel function. 
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Now we discuss briefly the dissipation effect on spin tunneling. For a spin tunneling 
problem, it is important to consider the discrete level structure. It was quantitatively shown 
that the phenomenon of MQC depends crucially on the width of the excited levels in the 
right well.S Including the effects of dissipation, the decay rate, in particular, is given byEffl 




where AE n is the level splitting, n' are the levels in the other well and tt nn > is the sum 
of the linewidths of the nth and n'th levels caused by the coupling of the system to the 
environment. For the exact resonance conditions, the temperature dependence of the decay 
rate is 



where the level broadening Q n contains all the details of the coupling between the magnet 
and its environment. If the width caused by the dissipative coupling sufficiently large, the 
levels overlap, so that the problem is more or less equivalent to the tunneling into the 
structureless continuum. In this case, the results obtained in this paper should be changed 
by including the dissipation. It is noted that the purpose of this paper is to study the spin 
tunneling at excited levels for single-domain FM particles in magnetic field at sufficiently low 
temperatures. Strong dissipation is hardly the case for single-domain magnetic particles,^ 
and thereby our results are expected to hold. It has been argued that the decay rate should 
oscillate on the applied magnetic field depending on the relative magnitude between the 
width and the level spacing.BlllfiBil However, it is not clear, to our knowledge, what 
should be the effect of finite temperature in the problem of spin tunneling. The full analysis 
of spin tunneling onto the precession levels remains an open problem. 





(28) 




(29) 



n 



B. tt/2 <9 h <tt 



For 7r/2 < djj < 7T, the critical angle 9 C is in the range of < 9 C < tt/2, and rj « 
y^2e/3. By applying the similar method, the problem can be mapped onto a problem of 
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one-dimensional motion by integrating out 0, and for this case the effective mass m and the 
effective potential U (5) in Eq. (f5) are found to be 

hS 2 (l + |cot^| 2/3 ) 



m 



2K 1 V 



l-e + 9K 2 (l + |oot 6 H \^) - K 2 '^f* + 2K 2 ^ e ^' 



(30a) 



and U (5) = 3U q 2 (l - jq), with q = 35/2^66, and 



U 



.3/2 _ 



I 4. /□ 1 1/3 
cot Vh\ 



2 7 ' 2 KiV 

3 3/2 h L 1 + jcot^H 



2/3 



15— 

1 K 2 

2 



(l + |cot^| 2/3 ) 



1/2 



(30b) 



The dependence of the effective potential E 1 (5) on 5 (= 9 — 9 ) for = 3ir/4 is plotted 
in Fig. 2, and (5) = 2K±VEi (S). The problem now becomes one of MQT, where the 
magnetization vector escapes from the metastable state at 5 = through the barrier by 
quantum tunneling. 

The nonvacuum bounce configuration with an energy E > is found to be 



tip — o v ^ [a - (a - 6) sn 2 (u^t, fc)] , 
o 



(31) 



where 



1 3C/ n 



^2 = 77 



a — c. 



2 V 2me 

a (E) > b (E) > c (E) denote three roots of the cubic equation 

E 



(32) 



o 3 n 

q ~ 2 q 



2U 



= 0. 



(33) 



sn^r, k) is the Jacobian elliptic sine function of modulus k = 

The classical action of the nonvacuum bounce configuration Eq. (31) is 

2 

+ U(S P ) =W + 2Ef3, 

J-/3 ^ \UT J 

with 







us p \ 


/ dr 




>-(3 







(34a) 
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w 



2 9/2 



^JmdT Q (a - cf 2 [2 (k 4 - k 2 + l) E (k) - (l - k 2 ) (2 - k 2 ) K (k)] . (34b) 



5 x 3 3 / 2 

The period t (E) of this case is found to be 

d5 



t(E) 



'2m 



2em 



-K (k') , 



(35) 



y/E-U (S) ^SU (a-c) 

where k' 2 = 1 — k 2 . Then the general formula Eq. (4) gives the imaginary parts of the 
metastable energy levels as ImE = exp (-W), where u (E) = ^|y, and W is shown in 
Eq. (34b). 

Here we discuss the low energy limit of the imaginary part of the metastable energy 

4 j> 



levels. For this case, E n = (n + \) fi 2 , Q 2 = \J ^U" {5 = 0) 



2V a ~ 2 1 1 



6 » (|A;' 2 )(l + fA;' 2 ) ) c 



A;' 2 (1 + iA;' 2 ), and A;' 4 = §§- < 1. After some cal- 



culations, we obtain the imaginary part of the low-lying metastable excited levels as 
hlmE n = (hlmEo), where 

2 25/4 x 3 5/4s e 5/4| cot ,1/6 



2 25/4 x 35/4^5/4 | CQt ^11/6 

<?2 — TJ7, ' 

1 6 4- ( 1 4- I not R l 2/3 ^ 4- 1 ft 7 2+9|cotg g | 2 / 3 3+|cot0 g | 2 / 3 

1-2 + 2^2 (^l + lcot^l j + 4^ (i+|cotej/|2/3)1 / 2 +^ 2(l+|coteH|2/ 3 ) 3/ 2 

The imaginary part of the metastable ground-state level is 

MmEo = (KlV)e '*s-v> ic ° tfe|i/ ; s [1 - \ + ^ (1 + Mfei") 

v 7 ^ 1 + |cot^| 2/3 l 4 4 v 1 ; 



4^4 

_1 51 -2|cot^| 2/3 1— 3 + |cot^| 2/3 
_ « 2 (l + |cotfe|^) 1/2 + 2 2 (l + |cot eH r) 3/2 

R pvnnnpnt it: 



exp (-W ) • 



(37a) 



where the WKB exponent is 



917/4 x 31/4 

W = Se 5/4 I cot 9 H I 

5 



1/6 



l + |-^ 2 (l + |cotfe|^) 1/2 



1 2 + 9|cot^| 2/3 — 3 + |cot^ 

'-^2 777T — -fV2 



4 (l+ |cot0 H | 2/3 ) 1/2 (l + |cot0# 

The decay rate at a finite temperature T is found to be 



2/3 



2/3 



3/2 



(37b) 



T (T) = 2Im£ (l - e - m2/3 ) exp (g 2 e~^ 2/3 ) 



(38) 
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In Fig. 3 we plot the temperature dependence of the tunneling rate for the typical values 
of parameters for nanometer-scale single-domain ferromagnets: S = 6000, e = 1 — H/H c = 
0.01, K 2 = 0.01, and 6h = Sn/i. From Fig. 3 one can easily see the crossover from 
purely quantum tunneling to thermally assisted quantum tunneling. The temperature Tjj ^ 
characterizing the crossover from quantum to thermal regimes can be estimated as ksT^ = 
AU/W , where AU is the barrier height, and W is the WKB exponent of the ground-state 
tunneling. It can be shown that in the cubic potential (q 2 — q 3 ), the usual second-order 
phase transition from the thermal to the quantum regimes occurs as the temperature is 
lowered. The second-order phase transition temperature is given by /cbTq 2 - 1 = where 
Ub = \J j^W" (xb)\ is the frequency of small oscillations near the bottom of the inverted 
potential —U (x), and Xb corresponds to the bottom of the inverted potential. For the 
present case, it is easy to obtain that 



k B T^ = - — {K X V) S-h 1 ^ 



2 l/4 x3 l/4 Ic0t^| 1/6 



^ " 1 • cot 0„ 2 3 



e 



9_ /„ , , „ l2/3 \V2 



i-- + -K 2 [i + \cote H 



1 21-2|cot^| 2/3 — 3+|cot^| 2/3 

:K 2 — ~~T7o + K 7 



4 + |cot0 H | 2/3 ) V2 (l + |cot^| 2/3 ) 3/2 

and k B T^ = (57r/18) ksT^ w 0.87/cbTq 2 \ For a nanometer-scale single-domain FM parti- 
cle, the typical values of parameters for the magnetic anisotropy coefficients are K\ = 10 8 
erg/cm 3 , and K 2 = 10 5 erg/cm 3 . The radius of the FM particle is about 12 nm and the 
sublattice spin is 10 6 . In Fig. 4, we plot the 9 H dependence of the crossover temperature 
T c for typical values of parameters for nanometer-scale ferromagnets at e = 0.001 in a wide 
range of angles 7r/2 < 9h < 7r. Fig. 4 shows that the maximal value of T c is about 0.26K at 
9 H = 1.76. The maximal value of T c as well as T is expected to be observed in experiment. 
If e = 0.001, we obtain that T c (135°) ^0.23K corresponding to the crossover from quantum 
to classical regime. Note that, even for e as small as 10~ 3 , the angle corresponding to an 
appreciable change of the orientation of the magnetization vector by quantum tunneling is 
5 2 = V6e rad> 4°. It is quite large enough to distinguish easily between the two states for 
experimental tests. 
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IV. MQC AND MQT FOR TETRAGONAL CRYSTAL SYMMETRY 



In this section, we study the FM particles with tetragonal crystal symmetry in a mag- 
netic field at arbitrarily directed angles in the ZX plane, which has the magnetocrystalline 
anisotropy energy 

E (9, 0) = K x sin 2 9 + K 2 sin 4 9 - K' 2 sin 4 9 cos (40) - M H X sin 9 cos - M H Z cos 9 + E , 

(39) 

where K 1: K 2 and K 2 are the magnetic anisotropy coefficients, and K x > 0. In the absence 
of magnetic field, the easy axes of this system are ±z for K 1 > 0. And the field is applied in 
the ZX plane as in the previous section. By using the dimensionless parameters defined in 
Eq. (8), and choosing K' 2 > 0, we find that = is an easy plane for this system, at which 
Eq. (38) reduces to 

1 



E (9, = 0) = | sin 2 9 + [K 2 - k'^ sin 4 9-Hcos(9- 9 H ) + E , 



(40) 



where K 2 = K' 2 /2Ki. Assuming that 
and the critical angle as 

1 



Ko-K, 



< 1, we obtain the critical magnetic field 



sin 9 r 



(sin9 H ) 2/3 + \cos9 H \ 2/3 
1 



3/2 



4 (k 2 - K 2 ) 

1 H — 

1 i l™+a__l 2 / 3 



1 + |cot0 H 



1/2 



1 2/3 



1 + 3(^-^ 2 ) 1 + i cot , ~m 



(l + |cot^| 2/3 ) 

Introducing 5 = 9 — 9 (\5\ <^ 1 in the small e limit), we derive the energy E (9, 0) as 



(41) 



E (5, 0) = K 2 [l- cos (40)] sin 4 (9 + S)+H x (l- cos 0) sin (0 O + 5) + #i (5) , (42) 



where E 1 (5) is a function of only S given by 
1- 



E 1 (5) 



-H c sin (9 C - Q H ) + (V 2 - sin (40 c ) 



(5 3 - 35 2 V ) 



+ 
+4 



cos (0 C 

o 



0tf)+ X 2 -^ 2 COs(40 c ) 



(5 4 - 45 3 ^ + 65V - 45 2 e) 



(K 2 - icQ e5 2 cos {49 c ) . 



(43) 
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A. H = tt/2 



For this case, we obtain that 77 rs \/2e 



1 - 4 (Kn - K. 



and 9 C = 7r/2. Performing the 



Gaussian integration over 0, we can map the spin system onto a problem of particle with 
effective mass m moving in the one- dimensional potential well U (5). For this case, 

hS 2 

m = 



2VK 1 



I- e + A[K 2 - K'\ +16i^ 



and 



U{5) 



Ah 



1 + 12 (K 2 - K, 



6 2 (6 - 2 V f 



(44) 



By applying the method similar to that in Sec. III. A, we obtain the low-lying tunnel 

splitting at degenerate excited levels as HAE n = (HAE ), where q 3 - 2 ~ ,s ' 

The energy shift of the ground state is 

2 13/4 



HAE = —=r {KxV) e 5/4 S~ 1/2 ( 1 - - A + AK 2 ) exp 
\/ir V 4 / 



-W a 



(45a) 



where the WKB exponent is 

2 5/2 



Se 3 / 2 



1 + 5 



8(K 2 -K 



8K n 



(45b) 



Eqs. (45a) and (45b) agree well with the result obtained with the help of the vacuum 
instanton solution.! And the final result of the decay rate at a finite temperature T is 
r (T) = (AE ) 2 [tt (1 - e' m ^) /2fi 3 ] I (2q 3 e' m ^/ 2 ) , where I (x) is the modified Bessel 
function. 



B. vr/2 < H < vr 



For 7r/2<0ff<7r, 77 ~ a/2c/3, the effective mass m is 

frS 2 (l + Icot^l 273 ' 



2#iV 



1 - e + 16if 9 + UK 2 - k'\ 3 - 2|coteH l 



,2/3- 



(46a) 



V l+|cot%| 2 / 3 J 
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and the effective potential is U (5) = 3U q 2 (l — with q = 35/2y/6e, and 



cot I 



HI 



1 1/3 



I cot ^| 2/3 



3 V V i + cot e H 2/3 



(46b) 



For this case, the imaginary part of the low-lying metastable excited levels is hlmE n 
% ( hlmEn), where o 4 = 2 / x3 / Se / |cot6> g | / — ^ imaginary part of the metastable 

ground-state level is 



5+8^2 + 51 ^2— ^2 I' — i T27T 



013/9 v 931/8 

Mm E = ^ (K.V) e^S- 1 / 2 

V 7r 



3Ot0 



1/4 



X 



1 + |cot^| 2/3 

12 loot 6 H \ 2/3 -7 



2 )~ 



+ I cot 6 H 



2/3 



exp (-W ) . 



(47a) 



where 



e — / 4 /_ _/\ |cot^| 2/3 -2 

1 + __ 8A - 2 __( A - 2 - A , 2 )L_£_ 



(47b) 



1 + \cot9 H \ 

The final result of the decay rate at a finite temperature T is T (T) = 
2Imi?o (l — e ~^ 4/3 ) ex P (?4 e_nn4/3 )- And the second-order phase transition temperature 
characterizing the crossover from quantum to thermal regimes is found to be 



(2) _2V4 x3 V4 /r _^_ ] i/4 |cot^| 1/6 



7T 



(K X V) S-h 1 ^ 



1 + |cot^| 2/3 



X 



e n — / 4/— — / \ 5 — 3 I cot 1 
l- 5 + 8A- 2 + -(AW 2 ) 7 



'if 



1 2/3 



+ I cot 6 H 



2/3 



V. MQC AND MQT FOR HEXAGONAL CRYSTAL SYMMETRY 

In this section, we study the hexagonal spin system whose magnetocrystalline anisotropy 
energy E a {6, 0) at zero magnetic field can be written as 



E a (9, 0) = K x sin 2 6 + K 2 sin 4 9 + K 3 sin 6 9 - K' 3 sin 6 9 cos (60) , 



(48) 



where K\, K 2 , K 3 , and K' 3 are the magnetic anisotropic coefficients. The easy axes are ±z 
for K\ > 0. When we apply an external magnetic field at an arbitrarily directed angle in 
the ZX plane, the total energy of this system is given by 
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E (9, 0) = E a (0, 0) - M H X sin 9 cos - M H Z cos 9 + £ , 



(49) 



By choosing K' 3 > 0, we take = to be the easy plane, at which the potential energy can 
be written in terms of the dimensionless parameters as 

E (9, = 0) = - sin 2 9 + K 2 sin 4 9 + (k 3 - k'^J sin 6 9 - H cos (9 - 9 H ) + E , (50) 



where K 3 = K 3 /2K 1 and K 3 = K' 3 /2K 1 . 

Under the assumption that \K 2 \, 
field H c and the critical angle 9 C as 



C 1, we obtain the dimensionless critical 



sin0 r 



(sin ^) 2/3 + | cos ^ 1 2/3 
1 



3/2 



1 + 



(l + |cot^| 2/3 ) 1/2 



i + 8_ loatfe 



1 + \cot9 H 

2/3 



2/3 



+ 



3 - K' 3 



(l + |cot^| 2/3 ) 2 



3 l + |cot0^ 



2/3 



+ 8 



s - K' 3 ) 



I cot 9 H 1 2/3 



(l + |cot^| 2/3 ) 

By introducing a small variable 5 = 9 — 9 (\5\ <^ 1 in the limit of e «C 1), the total energy 
becomes 



• (51) 



E (5, 0) = K 3 [1- cos (60)] sin 6 (9 + S)+H X (1- cos 0) sin (9 + 5) + E x (5) , (52) 
where E 1 (5) is a function of only 5 given by 

E x (5) = ^H c sin (9 C - 9 H ) + K 2 sin (40 c ) + 4 (k 3 - k'^J (5 sin 3 9 C cos 3 C - 3 sin 5 C cos 9 C ) 



x (5 3 - 35 2 r]) + pl c cos (0 C - 9 H ) + Z 2 cos (40 c ) + 3 3 - K' 3 ^j sin 2 C (sin 4 9 C 
-10 sin 2 C cos 2 9 C + 5 cos 4 C )] (<5 4 - 4<5 3 r/ + Q5 2 7] 2 - 45 2 e) + e<5 2 [4K 2 cos (40 c ) 
+12 (Z 3 -K'^j sin 2 c (sin 4 c - 10 sin 2 C cos 2 C + 5 cos 4 C ) . (53) 



A. H = tt/2 



For 0# = 7r/2, i.e., the external magnetic field is applied perpendicular to the anisotropy 



axis, we obtain that 9 C = tt/2 and r\ = v2e 



1 - AKo - 12 [K 3 - K 



. The spin system 
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can be mapped onto a particle with effective mass m moving in the one-dimensional potential 
well U (5), where 

HS 2 



m 



2VK 1 



l-e-AK 2 -6[Kn-KA- 3QK 



(54a) 



and 



U (5) = ^ [l + 12K 2 + 30 (K 3 - <)] 5 2 {5 - 2 V ) 



(54b) 



By applying the similar method, we obtain that the energy shift of the n-th excited level is 
HAE n = j (hAE ) ,where 

29/2^3/2 
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2 1 

The energy shift of the ground state is 

ol3/4 r / . , 

HAE = —j=- (KiV) e 5/i S- 1/2 |l - - - 6 \K 3 - + 9if 3 J exp (-W ) , (55a) 
and the WKB exponent is 

2 5/2 



I + - N?r, 24 - K'Q - i8/v; 



(55b) 



The decay rate at a finite temperature T is 



T (T) = (A£ ) 2 [tt (1 - e- m ^) /2ft 5 ] J fyse"* 1 ^ 2 ) , 



where 



Q 5 = 2 3/2^ e 3/2 

HS 



(k 3 - k' z ) + 18^;] 



l-- + 4K 2 + 6(ir 3 -^3) +18^31 



B. tt/2 <9 h <tt 



For this case, the effective mass m and the effective potential U (5) are 

hS 2 (l + |cot ^| 2/3 ) 



1 -e 



4"^ 3-2|cote g | 2 / 3 
oA 2 



i+|cote H | 



273" 



+ 2(K 3 -^ M 3_4|cot 



f |2/3 



i+|cote H | 2/3 )' 



36 ^3 1+ | cot% |2/3 
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and 



U(S) = 



1/3 



K X V \cot9 H 
ft l + |cot^| 2/3 

x5 2 (Vfk-S^ . 



1 + 4 7-4|catfa 



2/3 



- -,\ 11-16 |cot0 H 



3 l + |cot# H 



2/3 



+ 2 (jf 3 - 4) 



2/3 



(l + |cot^| 



2/3 



The imagi 



..aginary part of the metastable excited levels is hlmE n = ^(hlmE ), and the 

imaginary part of the ground state is 

» m£o = 37/9 «J 231/8 wys-v ' cot *"' 1/4 2/3 

v 7 ^ 1 + lcot^l 7 



' e 2 12|cot^| 2/3 - 
4 3 l + |cot^| 2/3 



2( __ ^9-13^*^ 



y - ia icotynr'" , CTF7 1 1 
(l + |cot^| 2 / 3 ) l + |cot^| 2 / 3 

where the WKB exponent is 



exp (—W ) > (56a) 



W = 2l7/ ' X3l/4 ^ 5 / 4 |cot^| 1 / 6 
5 



. e 4_ 2-|cot^| 2/3 

1 ~~ 7 + 7-^2 iJTT 

4 3 l+cot0 ff 2/3 



1 + |cot0 H | 2/3 
! — — ~ — 18-fTo ^ 

(i + |cote H | 2 ' 3 ) i + |cote H | 2/3 



(56b) 



and 



% = 



1 



c 4J7 2-|cot 6> g | 2/3 



2 25/4 x 35/4^5/4^0^1/6 
3 . /— — /\ 2~3|cotc9 g l 2/3 



4 TV- Z-JCOTPHJ ' A I TV" 7^' \ Z-JJCOtPHJ ' | 10^' i 

7 ' ' ' ~ ~^73 - 4 I tf 3 - * 3 J (l+ | cot ^/3) 2 + 18 ^3 



(56c) 



'l+lcote^ 3 -V" J "V (i+|cote ff | 2 / 3 ) z ' 3 i+|cote w | 2 / 3 
The final result of the decay rate at a finite temperature T is T (T) 
2ImE (l - e~ m(il3 ) exp (g 6 e _? ^ 6/3 ), where 



I V6 



fi 6 = 2 5 / 4 x 3 i/4^ e i/4^t^ 

^ l + |cot^| 7 



, e 4— 5-3|cot^ 
1 — - + -A 2 



2/3 



2 ' 3"' l + |cot^| 2/3 
1 



+ 2(K 3 -K' y- W ^Z +^ 

V ' (l + jcotfe| 2 ") 3 l + |cot9 H | 2 ' 3 

The second-order phase transition temperature characterizing the crossover from quantum 
to thermal regimes is found to be 
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k B T™ 



>l/4 



gl/4 



TT 



Icot 6»h 



1 1/6 



+2 



s - ^3 



1 + |cot^| 2/3 

7-10lcot e H \ 2/3 



e 4— 5-3|cot^| 2/3 
2 + 3 2 l + |cot^| 2/3 



(l + |cot^| 2/3 )' 



+ 18K» 



1 

'1 + Icot^ff 



2/3 



VI. CONCLUSIONS 

In summary, we have theoretically investigated the quantum tunneling of the magne- 
tization vector between excited levels in single-domain FM nanoparticles in the presence 
of an external magnetic field at arbitrary angle. We consider the FM particles with the 
general structure of magnetocrystalline anisotropy. By calculating the nonvacuum instan- 
ton in the spin-coherent-state path-integral representation, we obtain the analytic formulas 
for the tunnel splitting between degenerate excited levels and the imaginary parts of the 
metastable excited levels in the low barrier limit for the external magnetic field perpendic- 
ular to the easy axis (9 H = 7r/2), and for the field at an angle between the easy and hard 
axes (tt/2 < 9 H < tt). The temperature dependences of the decay rates are clearly shown 
for each case. The low-energy limit of our results agrees well with that of ground-state spin 
tunneling. One important conclusion is that the tunneling rate and the tunnel splitting at 
excited levels depend on the orientation of the external magnetic field distinctly. Even a 
small misalignment of the field with 9 H = n/2 orientation can completely change the results 
of the tunneling rates. Another interesting conclusion concerns the field strength depen- 
dence of the WKB exponent in the tunnel splitting or the tunneling rate. It is found that 
in a wide range of angles, the e (= 1 — H/H c ) dependence of the WKB exponent is given 
by e 5/4 (see Eq. (37b)), not e 3 / 2 for 9 H = tt/2 (see Eq. (21b)). As a result, we conclude 
that both the orientation and the strength of the external magnetic field are the controllable 
parameters for the experimental test of the phenomena of quantum tunneling and coherence 
of the magnetization vector between excited levels in single-domain FM nanoparticles at 
sufficiently low temperatures. If the experiment is to be performed, there are three control 
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parameters for comparison with theory: the angle of the external magnetic field the 
strength of the field in terms of e, and the temperature T. Furthermore, the 9jj dependence 
of the crossover temperature T c and the angle corresponding to the maximal value of T c are 
expected to be observed in further experiments. 

In order to avoid the complications due to distributions of particle size and shape, some 
groups have tried to study the temperature and field dependence of magnetization reversal 
of individual magnets. Recently, Wernsdorfer and co-workers have performed the switch- 
ing field measurements on individual ferrimagnetic and insulating BaFeCoTiO nanoparti- 
cles containing about 10 5 -10 6 spins at very low temperatures (0.1-6K).! They found that 
above 0.4K, the magnetization reversal of these particles is unambiguously described by 
the Neel-Brown theory of thermal activated rotation of the particle's moment over a well 
defined anisotropy energy barrier. Below 0.4K, strong deviations from this model are evi- 
denced which are quantitatively in agreement with the predictions of the MQT theory with- 
out dissipation.^ The BaFeCoTiO nanoparticles have a strong uniaxial magnetocrystalline 
anisotropy.! However, the theoretical results presented here may be useful for checking the 
general theory in a wide range of systems, with more general magnetic anisotropy. The 
experimental procedures on single-domain FM nanoparticles of Barium ferrite with uniaxial 
symmetry! may be applied to the systems with more general symmetries. Note that the 
inverse of the WKB exponent B~ x is the magnetic viscosity S at the quantum-tunneling- 
dominated regime T -C T c studied by magnetic relaxation measurement s.0 Therefore, the 
quantum tunneling of the magnetization should be checked at any 9jj by magnetic relax- 
ation measurements. Over the past years a lot of experimental and theoretical works were 
performed on the spin tunneling in molecular Mni2- Ac§ and Fe« clusters having a col- 
lective spin state S = 10 (in this paper S = 10 3 — 10 5 ). These measurements on molecular 
clusters with S — 10 suggest that quantum phenomena might be observed at larger system 
sizes with S ^> 1. Further experiments should focus on the level quantization of collective 
spin states of S = 10 2 -10 4 . 

The theoretical calculations performed in this paper can be extended to the AFM parti- 
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cles, where the relevant quantity is the excess spin due to the small noncompensation of two 
sublattices. Work along this line is still in progress. We hope that the theoretical results 
presented in this paper may stimulate more experiments whose aim is observing quantum 
tunneling and quantum coherence in nanometer-scale ferromagnets. 
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Figure Captions: 

Fig. 1 The S (= 9 - 9 ) dependence of the effective potential E x (5) for 9 H = ix/2 (MQC). 

Fig. 2 The 5 (= 9 - 9 ) dependence of the effective potential E 1 (5) for 9 H = 3tt/4 (MQT). 
Here, K 2 = 0.001. 

Fig. 3 The temperature dependence of the relative decay rate T (T) /T (T = OK) for FM 
particles in a magnetic field with a range of angles n/2 < 9 H < tt. Here, S = 6000, 
e = 1 - H/H c = 0.01, K 2 = 0.01, and 9 H = 3tt/4. 

Fig. 4 The 9h dependence of the crossover temperature T c for ir/2 < 9 H < vr. 
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